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Majorana fermions are rising as a promising key component in quantum computation. While the prevalent 
approach is to use a quadratic (i.e. non-interacting) Majorana Hamiltonian, when expressed in terms of Dirac 
fermions, generically the Hamiltonian involves interaction terms. Here we focus on the possible pair correlations 
in a simple model system. We study a model of Majorana fermions coupled to a boson mode and show that 
the anomalous correlator between different Majorana fermions, located at opposite ends of a topological wire, 
exhibits odd frequency behavior. It is stabilized when the coupling strength g is above a critical value g c . We 
use both, conventional diagrammatic theory and a functional integral approach, to derive the gap equation, the 
critical temperature, the gap function, the critical coupling, and a Ginzburg-Landau theory allowing to discuss a 
possible subleading admixture of even-frequency pairing. 


Introduction Recent advances in material fabrication 
techniques have brought Majorana fermions out of a theo¬ 
retician’s toy box and placed them onto the stage of realis¬ 
tic objects. Majorana fermions are particles with the peculiar 
property of being their own anti-particles. While its status 
as a fundamental particle is still unclear in the arena it was 
first proposed, namely neutrino in particle physics, in con¬ 
densed matter physics, on the other hand, there have been 
realistic proposals and indeed fairly convincing experimen¬ 
tal signatures that these will come into existence as zero en¬ 
ergy modes [1-8]. Apart from the obvious theoretical inter¬ 
est, such Majorana zero modes are particularly attractive as 
a potential key component in topological quantum computa¬ 
tion [9-11] because their origin is often of a topological na¬ 
ture (vortices, end modes, topological Dirac materials, etc.), 
hence their existence is robust and protected against disrup¬ 
tions such as defects, perturbations and decoherence. Since 
they also exhibit anyonic statistics under suitable conditions, 
their mutual braiding would constitute the fundamental oper¬ 
ation to encode and process information. They are thus well 
suited for fault-tolerant quantum computation. 

While Majorana fermions are most often studied using an 
effective Hamiltonian bilinear in Majorana operators, the fun¬ 
damental object in the condensed matter context, however, 
is still Dirac fermions. These effective Hamiltonians, when 
written in terms of Dirac fermions, generically contain both 
particle (c^c) and pairing ( cc ) terms, which is why any con¬ 
densed matter platform for Majorana fermions necessarily 
contains superconductivity [5]. As pairing inevitably arises 
from interaction, an important question one should ask is what 
types of instability would ensue. For conventional fermions 
and bosons, such instabilities would generate non-vanishing 
anomalous correlation functions which represent the onset of 
steady orders such as charge, spin, or pair susceptibilities. The 
same can be done for Majorana fermions [12, 13]. The effec¬ 
tive Hamiltonian formalism, being time-independent in its na¬ 
ture, implicitly assumes that these instabilities are dominated 
by their equal time behavior (i.e. the even-in-time compo¬ 


nent). However, the very nature of Majorana fermions be¬ 
ing simultaneously a particle and an anti-particle hints to the 
possibility that time (or frequency) dependence in the order 
parameter has a fundamental role to play. As a heuristic ex¬ 
ample, the pairing correlator M r = (Tfi(j)/i( 0)) of Majorana 
operator p separated by (Matsubara) time r is forced to be 
odd in time , and therefore vanishes at equal time. A more 
complete discussion on the pairing symmetries of Majorana 
fermions will be presented later. 

Historically, an odd frequency (odd-f) pairing state was first 
pointed out by Berezinskii [14] in 1974 as a candidate state for 
superfluid He3. While this was eventually proven not to be 
the case, Berezinskii’s pioneering paper set the stage for sub¬ 
sequent searches of other possible odd-f states, such as spin 
singlet odd-f pairing state [15, 16]. Although initial search of 
models and systems realizing such states turned out to be in¬ 
conclusive [17], more recent discussions of odd-f states have 
moved to systems of superconducting heterostructures where 
one can have conventional pairs converted into odd-f states at 
the interfaces [18-21]. Moreover, odd-f states can also be re¬ 
alized in homogeneous multiband superconductors [22, 23]. 
If we take the broader view that odd-f states represent a novel 
class of hidden order, then this concept may be extended be¬ 
yond superconductivity to encompass a whole set of states that 
might have odd-f order parameters [17-19], e.g., spin nemat¬ 
ics [24], BEC [25], and density waves [26]. 

In this work, we will study pairing of Majorana fermions in 
the presence of interaction induced by coupling to an external 
boson. Our focus will be on the pairing between different Ma¬ 
jorana modes because it arises solely as a consequence of in¬ 
teraction. Unlike same-mode pairing which, as discussed be¬ 
fore, is required to be odd-f due to fermionic statistics, there 
is no a priori requirement on the frequency dependence for 
cross-mode Majorana pairing. Indeed, as we will show, both 
even-f and odd-f solutions exist for the gap equation, yet the 
odd-f solution has a lower free energy and is hence more stable 
than the even-f state. We determine the phase diagram in the 
temperature-coupling constant plane showing the disordered 
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FIG. 1. Schematics of superconducting wires that host Majorana 
fermions p l a on their edges, where i is the wire index and a = 1,2 
labels the edge. Dashed line illustrates possible coupling between 
Majorana fermions on neighboring wires. 

high temperature phase and a pair-correlated low-temperature 
phase, separated by a second-order transition. Interestingly, 
we find a quantum phase transition at a critical coupling, in 
contrast to conventional BCS theory where any weak attrac¬ 
tive interaction induces pairing, if only at an exponentially 
small temperature. We derive the gap equation and solve it 
both, analytically in the limit of high and low temperature, 
and numerically. The critical coupling at low T may be linked 
to the first excited state of an effective Schrodinger equation, 
enabling a determination of the critical coupling. We also de¬ 
rive a functional integral representation of the partition func¬ 
tion in terms of the pair order field, which allows to estimate 
the condensation energy. The latter may also be used to derive 
an approximate Ginzburg-Landau expansion in terms of both, 
the odd- and even-frequency order parameters. 

Majorana pairing state and Extension of Berezinskii's clas¬ 
sification To set the stage, first consider the notion of pairing 
state. A state is paired when any two-particle anomalous ex¬ 
pectation is finite: for any particle operator O a (r,r ), where a 
denotes a collective quantum number, its pairing amplitude is 

P a £{r,T) = {TO a (r,T)0\ 0,0)) (1) 

The cases of fermions 0(r,r) = c(r,r), spins O a (r,r) = S a (r,r) 
(<a denoting directions), and boson O a (r,r ) = b a (r,r) («a denot¬ 
ing, say, species) all fall within the category of thus defined 
paired states. This does not necessarily mean there is a true 
superconducting or superfluid order. For example, a typical 
pairing state that supports a superflow would exhibit a f/(l) 
gauge symmetry breaking related to phase coherence and off- 
diagonal long range order, and has other attributes like phase 
stiffness. We instead accept that thus defined pairing correla¬ 
tion is the basic property that we investigate in the context 
of Majorana fermions. Since Majorana operators are real, 
the only symmetry left to be broken is Z 2 symmetry. Thus 
there are important distinctions between the “paired” Majo¬ 
rana states and conventional pairing. Yet the question about 
the appearance of anomalous order in models with Majorana 
states stands and needs to be addressed. We focus on the pair¬ 
ing correlator as a fundamental ingredient that any “super” 
state must possess - such correlations might be important in¬ 
dicators of incipient orders. 

We start with the extension of the Berezinskii classification 


for Majorana operators. Consider a set of Majorana states that 
emerge from materials like the ends of magnetic or semicon¬ 
ducting wires, see Fig. . We introduce Majorana operators 

A4 = a4) f . i=L..N,a=l,2, (2) 

where i is the wire index and a denotes the two ends of the 
wire. Following Eq. , consider the matrix of pair amplitudes 

< fo (r) = <774 (t)^(0)> = . (3) 

The antisymmetry follows from the definition of time order¬ 
ing, regardless of the Majorana nature of the operators and is 
similar to Berezinskii’s classification. Majorana fermion pair¬ 
ing is analogous to same spin Dirac fermion pairing. 

In the simplest case where i = j and in the absence of any 
interaction, one has 

Kb(r) = M ab (r) = M{r)6 ab = -M{-r)6 ab . (4) 

Thus the only pairing channel open to a single Majorana 
fermion is the odd-f pairing. Mixed pairs may form in the odd- 
f channel, with a possible admixture of even-f component. 
From the free action S = f driid T \i = 2 2 w „>o iu n Hw n H-co n , the 
propagator, which is also the pairing amplitude, reads 

= f— , M(t) = i-sgn(r). (5) 

2ioj n 4 

Thus a free Majorana state has odd-f pairing built in. Upon 
analytic continuation, M(oj + id) = M'(oj) + iM"{co ), with 

M'(co) = PfT , M"(oj) = ?-6(a >), (6) 

2 co 2 

where Pf stands for principal factor. Some comments are in 
order: i) this analysis suggests that a free Majorana state at 
zero energy is the simplest realization of odd-f state. Since 
the Majorana fermion is both a particle and a hole, the single 
particle propagator is a pair propagator. The natural question 
would be to identify the pairing states that would emerge from 
interacting fermions, which we shall address later. There are 
attempts to build in this direction [27, 28]. ii) Real part of 
the pairing susceptibility is odd in frequency and time while 
the imaginary part is even. This reversal of parities of real 
and imaginary parts with respect to the case of BCS (even fre¬ 
quency pairing) is expected [15]. The imaginary part M"(qj) 
that represents the decay of the pairs has a sharp peak at zero 
energy, which indicates the strong decay of odd-f pairing at 
zero frequency, i.e. at long times. The situation here is similar 
but not identical to the case of odd-f states found in supercon¬ 
ducting heterostructures where a peak in DOS and in M"(qj) 
is present and is associated with Andreev scattering [19]. 

While the analysis above is for zero-energy Majorana 
modes, the inherent connection between single particle propa¬ 
gator and pairing remains for the general case of a propagating 
Majorana branch of exictations, 

G(k, r) = M(k, r) = (Tp- k (r)p k m • (7) 

In the remainder of this paper, therefore, we will not concern 
ourselves with the effect of dispersion, but instead focus on 
how interaction affects the fate of odd-f Majorana pairing. 
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Model of interacting Majorana fermions The minimal 
model for investigating the effect of interaction is to consider a 
pair of Majorana fermions. We focus in the following on a sin¬ 
gle topological wire with particle-hole symmetric edge states. 
Let c, c ^ be the electron operators, then the corresponding Ma¬ 
jorana operators p a ,a = 1 , 2 , are defined as coherent superpo¬ 
sitions of particles and holes that satisfy 

C=p\+ip 2 , — /i\ — i/J-2 • (8) 


We note that {p a ir),pbir)} = \d a b- The electron pair opera¬ 
tor may be expressed as c(r)c(0) = pi(r)pi(0) - p 2 (r)p 2 ( 0 ) + 
ipi(r)p2(0) + /yU2(r)yUi(0). The point to note now is that 
whereas pair correlations of like Majorana fermions are cor¬ 
relations in time of a single Majorana fermion at either end of 
the wire, pair correlations of unlike fermions, i.e. a finite am¬ 
plitude of p\p 2 , describe nonlocal correlations between Majo¬ 
rana fermions at opposite ends of the wire. The latter ampli¬ 
tude is zero for free fermions, but will be generated by inter¬ 
action, as we will show below. 

We now couple the Majorana fermions to a single boson 
mode of frequency Q, as described by the Hamiltonian 

H e b = gi^+b)i(jiui2-H2Vi) + Q(b^b+lj) , (9) 


where g is the coupling strength. Such a boson mode can 
arise from an optical phonon, for example. We have cho¬ 
sen to couple the boson equally to both particle and hole, as 
i(p\p 2 -yu 2 jUi) = c^c- \ - \c*c — \ccK The Euclidean ac¬ 
tion corresponding to Eq. 9 is S e b = f drf/d T il'r + f(d T + + 

gff + f)fnl /, where if/ and f are fields for c and b , respec¬ 
tively. Integrating out the boson then leads to an interact¬ 
ing action of the fermion, S = -\ogJr ( f ) e~ Seb = f drijfd T if/- 
g 2 f djdr' fjj\fj T K T yfj T of T ^ where K TT > - \{T(b + b^)(r)(b + 
^)(r')) is the boson propagator in time domain, 



-nir'-ri 2cosh[D(r'-T)] 
e e^-\ 


( 10 ) 


and (3 = T 1 is the inverse temperature. In the Majorana basis, 
using the notation p = p\,q = p 2 , the effective action is 


-/ 


dr(pd T p + qd T q) - £ 


1 J* dr dr' 


HtUt' drdr' • ( 11 ) 


Since K T y = K r > jT , the r > r' interaction is identical to that of 
r < r' and is accounted for by doubling the prefactor. 

Gap equation for pq pairing We now sketch out deriva¬ 
tion of the gap equation. To decouple the four-Majorana in¬ 
teraction, we employ Hubbard Stratonovich (HS) transforma¬ 
tion, after which the interaction term becomes (cf. SM) 


o HS 
0 int 


/ 


drdr' 


A t? t'At',T 

16 




( 12 ) 


where A r y are complex HS fields subject to the constraint 


A t?t ' 




(13) 


which originates from the aforementioned redundancy of the 
r > r' interaction (cf. SM), and the frequency domain con¬ 
straint follows assuming A depends only on r' - r. [29] Inte¬ 
grating out p and 77 , one obtains the effective action of A, 


co m ,co n 6 r 


-Z log 

Wm-Wn (jj n >0 


Cdi- 




(14) 


where o) n = (2 n + 1 ) 7 tT are fermionic Matsubara frequencies. 
Demanding SS Eff /dA^ = 0 then yields the gap equation, 


This equation may be expressed in the more familiar form 

A„„=8 (16) 

% 

where M 21 = (A - - A + )/[4 co 2 + (A - ) 2 - (A + ) 2 ] is the anoma¬ 
lous Green’s function (cf. Eq.3) found by inverting the matrix 

M~ l = 2iu) n - cr x A^ n - i(r y A^ n . (17) 

where A% n = \{A Un ±A _ w „), and the phonon propagator K Vn = 
Q/[(Q) 2 - ( iv n ) 2 ] is the Fourier transform of Eq. 10. Here 
we used the symmetry property of Eq. 3. It is worth men¬ 
tioning that the gap equation Eq. 16 differs from the conven¬ 
tional gap equation by a factor of 4. This is due to the fact 
that the diagram rules for Majorana fermions as compared to 
usual fermions have to be modified by applying a factor of n at 
each interaction vertex involving n Majoranas (at the electron- 
phonon vertex this would be a factor of 2), see also SM. 

Emergence of odd-f pairing We now discuss the solutions 
to the gap equation. The normal state (A^ n = 0) is a trivial 
solution. Nonzero solutions can be found analytically at high 
and low temperatures. At high temperatures T » Q, the ker¬ 
nel may be taken to be diagonal, = d mn /£l = d mn K T , T /T 

where K TT is the equal time boson propagator, see Eq. 10. 
Eq. 16 only has pure odd-f and even-f solutions, 

= sA ( ll n = 2 fl,flfr + sojI , s = ±1 . (18) 

We emphasize that no assumption was made on how A trans¬ 
forms under cj n <-> -o) n , yet the saddle point solutions acquire 
a symmetry. The pure odd-f (even-f) solution is found to have 
a lower (higher) free energy than the unpaired state (cf. SM), 
hence we will focus on the odd-f case. Stability of the odd-f 
state is consistent with Ref. [30] and is discussed in more de¬ 
tail in SM. The reality constraint Eq. 13 imposes a frequency 
cutoff oj c = sjlg 2 T /Q for the odd-f solution above which the 
gap function falls back to zero. Since the smallest possible 
\cj n \ is coo = nT , requiring o) c > oo entails that there is a criti¬ 
cal temperature T c above which the system is normal, 



A-,, 


J GJn —(L>m 


IA A 

4 t±0) n ^-OJn 


(15) 


■a>„ 



A w „ 6 R, 


(19) 
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FIG. 2. (Color online) Pairing of different Majorana fermions. Cou¬ 
pling strength is set to g = 6Q. The gap function peaks at the larger 
of the two energy scales, temperature nT or boson frequency Q. 



FIG. 3. (Color online) Phase diagram in log scale at Q = 1. Colored 
crosses: numerical results for the phase boundary T c (g). Green solid 
line: high T limit Eq. 19. Gray dashed line: T —> 0 limit Eq. 22. 


Closely below T c , the first nonzero gap component A Wo has 
the usual mean field temperature dependence 

A^=27rT c y/l-T/T c . (20) 


At low temperatures T «D, the kernel K Vn is approxi¬ 
mately constant for v n < Q and drops as v“ 2 beyond. Therefore 
we expect A^ n to be maximum at oj n ~ Q. Approximating A 
by its value at this typical frequency on both sides of Eq. 16, 
one is left with the frequency summation T Yv n K Vn = K TT ~ \. 
The resulting A^ n= ^ is again given by Eq. 18, hence a real val¬ 
ued solution can only be found for g > g c ~ Q. 

For generic temperatures, we have obtained the gap func¬ 
tion numerically by minimizing the action Eq. . Fig. 2 
shows results of ^ versus ^ for a fixed value of ^ = 6 and 
several values of T <T C . The peak sits at the dominant energy 
scale: tu pea k = nT for T » Q, and tu pea k ~ Cl for T —> 0. The 

slope ^ serves as an order parameter for odd-f pairing. 

Critical coupling and effective Schrodinger equation The 
T Cl analysis above indicates the possibility of a critical 
coupling g c ~ Cl at zero temperature. This agrees with our 
numerically determined phase diagram Fig. 3 . To clarify this 
question, we now analyze the gap equation at T = 0 and A —> 0, 
which may then be linearized and transformed into a differen¬ 
tial equation in time 


C v(x) 


Y(x) = 0 . 


( 21 ) 


where x-t'-t , Y(x) = —Y{—x) is the Fourier transform of 
A“/tu 2 , and V(x) = g 2 e~ a ^ x K This is a “Schrodinger equation” 
for energy eigenvalue E — 0. The condition for the lowest- 
lying odd-parity state to have energy E = 0“ is found as 


g c (T = 0) = -Q , (22) 

r 

where y - 2 if V(x) is approximated by a square well of depth 
g 2 and width Q -1 , in which case Y(x) oc sin(g c x) inside the 
potential well and decays rapidly on the outside, hence A~ n 
peaks at tu pea k = g c - Numerically we found y ^ 2.6, see Fig. 3. 
We have a quantum critical point at g = g c separating the nor¬ 
mal phase from the paired phase. 


Ginzburg-Landau theory In order to understand the inter¬ 
play of odd-f and even-f pairing better we now formulate a 
Ginzburg-Fandau free energy expression. As a suitable order 
parameter for odd-f pairing we choose the slope A of the fre¬ 
quency dependent gap function, A - = Aoj in the low frequency 
regime \oj n \ < co x . The even-f order parameter we define as 
B = A + /Q. Then, the free energy of low frequency modes as 
derived from Eq.14 may be approximated as (see SM) 

A 4 

F = c 0 [(f ■- 1) A 2 + — ] + [(bot + 4ci)B 2 + 2C2-B 4 ] - CiA 2 B 2 

(23) 

where t = b { = T c Z' m , n 16 ( ^-i ^ ] mn , \ 

and YJ is restricted to \aj n \ < oj x , with cj x chosen such that 
b\ = co- It is seen that as found numerically odd-f pairing 
appears for T <T C (provided g > g c ), whereas even-f pairing 
always increases the free energy, even taking into account the 
negative contribution of the mixed term (for details, see SM). 

Conclusion We extended Berezinskii’s classification of 
pairing order parameters to Majorana fermions, and showed 
that odd frequency pairing arises naturally in both free and in¬ 
teracting Majorana theories. For free Majoranas, odd-f pairing 
is the only pairing channel allowed by the fermionic exchange 
statistics. With interaction, pairing of Majorana fermions at 
opposite ends of a topological wire may be induced. We show 
that within a model of phonon-mediated interaction there is a 
critical coupling strength for pairing, which implies the ex¬ 
istence of a quantum critical point, and can be determined 
by mapping the gap equation onto an effective Schrodinger 
equation. We find possible pairing solutions both analytically 
and numerically and show that only odd-f pairing is stable, al¬ 
though even-f pairing is in principle allowed. To this end we 
define two frequency-independent order parameters character¬ 
izing both odd-f and even-f pairing and propose an effective 
Ginzburg-Fandau free energy derived from the microscopic 
mean field action. The simple model study of Majorana pair¬ 
ing presented here opens the way for a study of more complex 
systems involving many Majorana states, e.g. located at oppo¬ 
site edges of a two-dimensional topological superconductor. 
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Note added We recently became aware of a related work 
by Zhang and Nori [31]. 
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Supplementary Material 

In this note, we derive the gap equation in detail using both field theory and diagrammatic expansion. Then we discuss the 
free energy of the odd-f and even-f pairing fields. 


Field theoretic derivation of the gap equation 


The Majorana action from the text is 

S = So + Sint > So = J" dr(jdd T /j + rjd T rj) , 5mt = -8g 2 Jdrdr'K T y 


li T li T ’r] T T} T ' 


Kry = -\e 




2 cosh [Q(r' -r)] 


eP Q - 1 


K v = 




Q 2 + v 2 


(A.l) 

(A.2) 


Since K r y = K T ^ T , the r > r' interaction is equivalent to that of r < r' and is accounted for by doubling the prefactor in S\ n \. 
The purpose of keeping track of such redundancy in the interacting action is to avoid introducing multiple Hubbard-Stratonovich 
(HS) fields for equivalent interaction terms. 

HS transformation uses the following identity: for a real positive number M and any four Grassmann variables 01,03,03,01, 


e Mif/nf/ 2 if/3^4 - x 


f 


dAdA* e~^ A *~ Mlf/l ^ M 1( A -MfoM 


■f 


dAdA* exp 


M 


+ A^i^ 2 + AV3^4 


(A.3) 
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upon proper normalization of the integration measure dAdA*. To decouple S\ n \ in the jdrj pairing channel, we use Eq. A.3 and 
introduce, for each pair of time indices t,t' with r < r', a complex field A T y. After HS transformation, S j n t becomes 


,HS 
> int 


= J* drdr' 


|a tj 


8 g 2 K Ta 


■ A T T 'ii T rj T ' - A* T ,/v ^77- 


Note that A is only defined for r < r'. One can introduce a new field A r?r ' with unrestricted time indices, 

I A r?r / t < r' 


A t?t ' - • 


A*, T > T r 

, r ,r 


If A t?t ' only depends on the time difference r' - r, then Eq. A.5 entails that the Fourier transform of A(r' - r) is real, 

A^ n = J*drA(r)e l0JnT = J* dr(A(r)e l(x)nT + c.c.) < 

Using the A fields, Eq. A.^ can be rewritten as 


T> 0 


,HS 
1 int ' 


16,r^ r ,r' 

Reinstating 5 o and going to frequency domain, the full action is 


J* drdr' 


-A T yii T ri T f 


=G~ 


i HS 


2 


A&J m ^C0 n 


\6g 2 p 


+ J](Mco n , w a *_ 6 


°n-Um oj n > 0 


2i(l) n A^ \ If~l-CL) n 
■oj n 2iaj n ) yl-Qj n 


=S Fer 


The fermion part is bilinear and can be integrated out, 

Z Fer = f D(^)e~ S ^ = P| 


O) n >0 


--Det G~' 
4 


- P[ ( w n . A W( 1 A_ W(1 ) , 


io n > 0 


(A.4) 


(A.5) 


(A. 6 ) 


(A.7) 


(A. 8 ) 


(A.9) 


where the factor -1/4 comes from the Jacobian d(nt 0 ,id-t 0 ,rit 0 ,r]- w )/d(i{s a „ifr a „if/- (0 ,i{s-t 0 ). The effective action of A thus reads 


S E ff = E 


A<T> m ^C0 n 


16 g 2 P 


~Yj ic, g( w « - 


COm-tOn oj n >0 


and minimization with respect to A yields the gap equation 




A/,i — 


"CJ m 1 


IA A _, ,2 

69 m L "J CJ n -CO m 4 LS ‘O n LS —(x) n U 


(A. 10) 


(A.ll) 


Gap equation from diagrammatic expansion 

In this section, we treat the interacting action S\ n \ (Eq. A.l) as a perturbation, and derive the gap equation as the one-loop 
correction to the self energy. We will use v and oj to denote bosonic and fermionic Matsubara frequencies, respectively. 

Let us introduce two-component Grassmann vectors, 


X *'co = 


Vao 

dcO 


9 ^CO — 9 d~Co) 


(A. 12) 


The matrix Green’s function is defined as 


G<, = CPA) = (rp~" pp"). 

'\d<Jj/^-<JJ dcod-CO ) 


(A.13) 
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It is equivalent to the Majorana pairing matrix M defined in the text. One has 

*afi _ <'P ^- Sint >0 

where <• • • )o denotes averaging over the free action So, and in the vector notation, the interacting action can be written as 
-Sint =g 2 JcItcIt'K tt ' [^ r >F r | pV/iTV] 

= g 2 T K * PV ft?**] [^3^1 Kan -0)2- v)6(a>. 3 ~m + v), 

v."l, 2,3,4 


(A. 14) 

(A.15) 
(A. 16) 


where the matrix X = -cr y . The summand corresponds to the diagram Fig. 4a. To the order of g 2 , Eq. A. becomes 


Gl P = ^ ^ inl)) ° = <n^>0 + (XXi X (-S int)>o + 0(g 4 ) , 
I — W int/0 


(A. 17) 


where the superscript c denotes connected contractions, i.e. both and are to be contracted with fields in S j n t, not with 
each other. Explicitly written out, one has (summing over repeated indices) 

<'F“'i^x(-5int))o =g 2 T 2 X 6 (o>i-o) 2 -v) 6 (m-o> 4 +v)X a n cd x(VZxKy„ 2 % 3 X„ 4 x¥j c 0 . (A.18) 

V.W1.2,3,4 

If we ignore the Hartree terms in (• • • )£, i.e., those with the two external fields v F a and contracting on the same A vertex, then 
it corresponds to the diagram Fig. 4b, and since there is no distinction between particle and hole for Majorana fermions, can 
contract with any of a,b,c,d, and X F S can then contract with two of a,b,c,d on the opposite vertex of v F a . There are 4 x 2 such 
contractions of identical results. For example, consider the contraction of (a,b)(a,c)(d,/3). The two “anomalous” ones are 


/ \jja\jjb \ M0)ab /\t ja me \ r M0)ac 

2^(0 > T gj 3 /0 — • 

From the delta functions one has cj 4 = - 0 J 2 = oj and 0 J 3 = ~( 0 \ =oj-v. The contribution to the RHS of Eq. A. 18 is thus 


l T Y Kyl ab l cd x =g 2 T 


G^l\Y K ' a ' 


( 0 ) 

oj-v 


AG 


( 0 ) 


af3 


(A.19) 


(A.20) 


where the minus sign comes from exchanging the order of a and b in the contraction, and we used -A r = A. Eq. A.18 is then 


where 


<n^x(-s lnt )>g= [g^g^ 

£ w = 8xg 2 r2^vpGLV] • 


(A. 21 ) 


(A.22) 


Comparing Eq. A. 17 with the definition of self energy, 


^a> — Ww ) 


■Gd G,, 


: Gf +Gi? ) E, 


/^( 0 ) , /^( 0 )y /^( 0 )y ✓AO) , 


(A.23) 


one can identify that up to leading correction, Eq. A.22 is indeed the self energy. Replacing G^ with the full G on the RHS of 
Eq. A.22 then leads to the self consistent equation for the self energy matrix, 

= 8 X g 2 T Y Ky [aG oj -vA\ • (A.24) 

y 


Now writing 


A-i = 1 2 ico -A 0J \ 
" -\A-oj 2 ico) 


Go = 


/ 2 ico A w \ 
\-A_ w 2 ico) 

DetG - 1 


AGojA = 


2 ioj 

A/,i 


A-^ 

2 ico 


DetG; 


-1 


(A.25) 









(a) Interaction vertex 


(b) Self energy 


FIG. 4. Feynman diagrams, (a): Interaction vertex, see Eq. A. 16. Frequencies going into the same A vertex sum to zero, and arrows denote the 
sign of the corresponding frequency in the 6 function, (b): Self energy at the one loop level, see discussion below Eq. A. 18 for multiplicity. 
Thick line represents the dressed Green’s function used in the self consistent equation Eq. A.24. 


note that G 1 is the same fermion kernel as one would get in the field theoretic approach, cf. Eq. A. 8. Then the self consistent 
equation for the matrix element reads 


a w = 8 


A„ 


W-CjJ' 


A^A-a/ -4 co' 2 


(A.26) 


which is the same as Eq. A.l obtained from minimizing the HS action. Eqs. A.l and A.26 appear to have different factors of 
T, this is due to the convention of the Fourier transform A v = f drA r e lVT , and the fact that K r is dimensionless, thus both K v and 
(1 /K) v have a dimension of time, hence the different factors of T. 


Free energy consideration of the odd-f and even-f pairing states 


Consider the action Eq. A. 10 relative to the normal state A = 0 (which is always a solution of the gap equation), 


<55 (A) ee 5 E ff (A) - 5 Eff(0) - V 


16 g 2 /3 

(o m ,aj n ° r 


- 2 , lo gi 1 

OJm-OJn QJ n >0 



(A.27) 


For any even-f configuration A c , n = A a , n , 6S (A) > 0, because the first term is positive definite, and in the second term, log(l - 
A2 f! /4m2)<0. Thus without solving the gap equation, one can already conclude that any even-f state will be unstable. 

In the diagonal limit K Vn = S Vn ^Ko, the odd-f solution to Eq. A.26 is 


= -AI<u„ = 2 y[g^ Ko-a%. 

From Eq. A.27, the action contributed by and A_ W/? is 


OJn 


2 g 2 TK 0 


-log 


2jTKo 


(A.28) 


(A.29) 


One can verify that 1 -x + logx < 0, with equality occuring at x = 1, thus the odd-f solution Eq. A.28 has a lower free energy 
than the normal state and is stable. 

To understand in a more general way why only the odd-f state is stable, one may adopt the Ginzburg-Landau perspective and 
expand the log term in Eq. A.27 to quartic order in A, keeping both, odd-f and even-f components, A% n = ^(A^ n ± A_^ n ), 


SS( A) 


+ ^qj m ^qj n 


^ 16 g 2 p 

Um,U n ° ^ 


* 1 . 

Wm-Wn (jj n > 0 


(Al) 2 - (A- ) 2 (Al) 4 + (A- ) 4 - 2(A+ ) 2 (A- ) 2 


4 col 


32 m 4 


+ 0(A 6 ). (A.30) 


The first term on the r.h.s. does not have a mixed term owing to the even symmetry of 1 /K. We now immediately see that the odd- 
f component contributes a negative contribution in quadratic order, which may induce a transition, while the even-f contribution 
does not. However, while the fourth order terms of both pure odd-f and pure even-f type are positive, the mixed term is negative. 
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This opens the possibility that in an odd-f paired state, even-f pairing may be induced at some lower temperature. In order to 
explore this scenario we first simplify the analysis by defining frequency independent dimensionless order parameters A,B , 
where A is the initial slope of A - , A~ n = Aoj n and B = A^ /Q . We restrict the frequency integrations to low frequency, \co n \ < co x , 
where the cutoff oj x is to be suitably defined later. The above action may then be expressed as a Ginzburg-Landau free energy 
expansion in the order parameters A, B , 


F(A, B) = SS (A) = c 0 


( t-l)A 2 +— 


+ [( b 0 t + 4ci) B 2 + 2C2B 4 ] -c x A 2 B 2 + 0(A 6 ,B 6 ,---). 


(A.31) 


Here we defined t= jr and the positive dimensionless coefficients 


b t = T c J] 






i6g 2 n 2(f -'> 


Cl 


=-x 

0J n > 0 


Q 


2 <jJ n 


It 


(A.32) 


the primed sum imposes the aforementioned frequency cutoff oj x implicitly determined by requiring 


b\ = c 0 


(A.33) 


If we now substitute the meanfield solution A 2 = 4(1 - T/T c ) for odd-f pairing (assuming g > g c ) into the mixed term we find 
that the determining equation for B, 


—^ = (bo + 4c l )t + 4c 2 B 2 = 0, 
d(B z ) 


(A.34) 


does not have a real solution for any positive t. 



